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Abstract 

We study the semi-inclusive hadron production in deep inelastic scattering at small x. A trans- 
verse momentum dependent factorization is found consistent with the results calculated in the 
small-x approaches, such as the color-dipole framework and the color glass condensate, in the 
appropriate kinematic region at the lowest order. The transverse momentum dependent quark dis- 
tribution can be studied in this process as a probe for the small-x saturation physics. Especially, 
the ratio of quark distributions as a function of transverse momentum at different x demonstrates 
strong dependence on the saturation scale. The dependence of the same ratio is also studied 
by applying the Collins-Soper-Sterman resummation method. 
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There have been compeUing theoretical arguments and experimental evidence that satu- 
ration M plays a very important role in high-energy hadron.c scattering processes, and a. 
effective theory called Color- Glass- Condensate emerges to describe the relevant physics [3!,|4|. 
In particular, the parton distributions at small-x (x is the longitudinal momentum fraction of 
the hadron carried by the parton) and/or of large nucleus can be calculated from this effective 
theory, and they all demonstrate a saturation behavior. The rapidity {Y = Inl/x) evolution 
of these distributions are controlled by a nonlinear JIMWLK equation 
been thoroughly studied in the last decade. By employing the saturation physics, the deep 
inelastic scattering QDIS) structure function measured by the HERA experiments can be 
very well described js], y, Q, as well as the diffractive structure functions [ll, 12, Q 
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Forward hadron suppression in ci -|- Au colli- 



sions at RHIC experiments also indicates the importance of saturation in the small-x re- 
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gion [18|, Il9|, l20|, l21|. All these successes have encouraged rapid developments in small-x 



physics in the last few years [22 1. 

One of the key predictions of this effective theory is the transverse momentum depen- 
dence of the parton distributions in big nucleus at small-x, especially the gluon distribu- 

nnn 

tion [3|, l23|, l2J] . In the inclusive DIS process, the gluon distribution is convoluted into a dipole 
cross section, which only provides an indirect probe. In this paper, we argue that the trans- 
verse momentum dependent parton distributions can be directly probed in semi-inc 



processes, for example, in the semi-inclusive hadron production in DIS (SIDIS) 251]. In 



usive 



this process, there are separate momentum scales: the momentum transfered by the virtual 
photon squared and the transverse momentum of the observed hadron in the final state 
p±. Because of the additional hard momentum scale Q^, the final state hadron transverse 



momentum can be directly related to that of the parton distribution in the nuc 
when is much larger than p^. The relevant QCD factorization theorem 26 



eon/ nucleus 
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28| has 



been rigorously studied for the leading-power contribution to the differential cross section. 
In the following calculations, we will extend this factorization argument to the case that 
involves saturation physics, and we argue that the transverse momentum dependent factor- 
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ization formula is still valid in the so-called geometric scaling regime 

is much larger than the saturation scale Ql, but saturation effects are still important. 
As an example, we will demonstrate this factorization for the semi-inclusive DIS process at 
small-x. 
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FIG. 1: Semi-inclusive DIS at small-x, where the cross represents the quark fragmenting into the 
final-state hadron. The quark carries the momentum fraction ^ of the virtual photon, and transverse 
momentum k±. 

In the semi-inclusive DIS process, 

e + piA) ~^e' + h + X , (1) 

we observe the final state hadron with characteristic kinematic variables, such as the lon- 
gitudinal momentum fraction of the virtual photon and transverse momentum p±. The 
usual DIS kinematics variables are defined as = —q.q, xb = Q"^ /'^Pa ■ 1, y = I' PaIC-- Pa, 
and Zh = Ph ■ Pa/q ■ Pa, where Ph, i, Pa and q are momenta for the final state hadron, in- 
coming lepton and nucleon (nucleus), and the exchanged virtual photon, respectively. The 
transverse momentum p± is usually defined in the center of mass frame of the virtual photon 

ig. 1, we show the schematic diagram for this process in the 



and the incoming hadron. In 
dipole framework at small- a; 2^, where the virtual photon splits into a quark-antiquark 
dipole, then scatters off the nucleon/nucleus target before the quark (antiquark) fragments 
into a final state hadron. In the current fragmentation region (forward direction of the 
virtual photon), the quark-fragmentation contribution will dominate the cross section. 
The differential cross section for the above process can be calculated in the dipole for- 



malism 
have 



24| or in the classical Yang-Mills effective theory approach 23|], and we readily 



da{ep ^ e'hX) oil^Nc sr^ 2 f dz D{z 



where D{z) is the quark fragmentation function into the final state hadron, F{q±,XB) the 
unintegrated gluon distribution defined below, = Zh/z, and the fragmenting quark's trans- 
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verse momentum is k± = p±/z} The variable b here is a suppressed variable in F{q±,XB) 
which is defined as the impact parameter with respect to the center of the nucleus. If 
one assumes that the nucleus is cylinder like and nucleons are uniformly distributed inside, 
one can easily see that the b dependence is trivial and evaluate J dPb which yields vri?^ 
with R being the effective radius of the nucleus. The phase space factor dV is defined as 
dV = dx sdQ"^ dzhdp\, and Ti reads as 

2 



ki. k± - q± 



+(i-.)4f(i-o-Q-(^- (,^_;^).^^. ) , (3) 

where = .^(1 — CjQ"^- We have also taken the massless-quark limit in the above formula for 
simplicity, and the first term is the contribution from transversely polarized photons while 
the second one corresponds to longitudinally polarized photons. The unintegrated gluon 
distribution is defined through the Fourier transform of the dipole cross section: 

where Tqq is the scattering amplitude, and is characterized by the saturation scale which 
depends on x. This unintegrated gluon distribution contains the saturation physics, which 
diagrammatically represents the multiple scattering of the quark- ant iquark dipole on nu- 
cleon/nucleus target. When integrating over transverse momentum p± and the fragmenta- 
tion function using Ylih i dzDq^h{z) = 1, the above formula will reproduce the total DIS 
cross section in ep{A) — * eX ^. 

In this paper, we are interested in the factorization property of the above differential 
cross section in the kinematic region where is much larger than the final-state hadron 
transverse momentum p\. In the current fragmentation region, Zh is of order 1. Therefore 
the quark transverse momentum k± is of the same order as p±. Furthermore, we assume 



^ We could add transverse dependence to the fragmentation function D{z) and write it as D(z,pi±). In 
this case, 'H{^,k±) should also have pi± dependence and becomes 7i(^, P-^~p^^ ) accordingly. Here we 
approximately neglect the transverse dependence of the fragmentation functions. The transverse depen- 
dence may change the numeric results, however, it will not affect the following factorization discussion 
once the above change of variables is made. 

^ In the fcf-factorization at small- a;, the gluon momentum fraction x differs from xb because of the kinematic 
constraints ,33.] . In the leading logarithmic (Inl/x) approximation, these two are consistent. 
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that is also much larger than the saturation scale Q\ which sets the transverse momen- 
tum of the unintegrated gluon distribution. Under these limits, we will be able to study 
the transverse-momentum-dependent factorization, where we can separate the transverse 
momentum dependence of the final-state hadron into the incoming quark distribution and 
fragmentation function and/or soft factor 26|, |28|. An important advantage to utilize the 
above limits is that we can apply the power counting to analyze the leading power contri- 
bution, and neglect the higher order corrections in terms of p^/Q"^ where p± stands for the 
typical transverse momentum ~ A;_|_ ~ g^). 

Moreover, we notice that the integral of Eq. (2) is dominated by the end point contribution 



of ^ ~ 1 where is in order of k\ 2J] . In order to extract the leading power term from this 
equation, we can introduce a delta function in Eq. (2): J d^5{S, — 1/(1 + A^/e^)) ^, where 
A^ = (1 — + 0^-L ~ Q^Yi integrate out ^ first. This delta function can be further 
expanded in the limit of p]_^ Q^, 
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A2 
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(5) 



i \ (1-0 (1-0 

where a logarithmic term in the above expansion is power suppressed and has been neglected. 
The contribution from the second term is also power suppressed. To see this more clearly, we 
can substitute = ^A^/ (1-0 into Eq. (3), and the hard coefficient Ti will have an overall 
factor (1 — 0^- Combining this with the delta function expansion, we will find that the 
second term is the above expansion is power suppressed relative to the first one. Applying 
the delta function expansion in Eqs. (2) and (3), we will obtain the leading contribution to 
the differential cross section in the limit of p±_ <^ Q, 



da{ep — >■ e'hX) 



2vr3g4 Z^^f 



Xb 



X 



dVq^F{q_L, XB)A(g_L, k±) 



where 



A{q^,k^ 



kL\k± - q±\ 



k± - q± 



{l-^)kl + ak±-q±y \k^-q±\ 



(6) 



(7) 



^ If we replace the gluon momentum fraction xb by a; = xb/^ in Eq.(2), we will reproduce the kt- 
factorization formula 



33| with this delta function. 
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We noticed that the longitudinal photon contribution is power suppressed and has been 
dropped. 

On the other hand, a transverse momentum dependent factorization can also be used 
to describe the SIDIS process when the hard scale (Q^) is much larger than the trans- 



verse momentum scale p\. To leading power of p^/Q"^, for example, we will have 



factorization formula for the differential cross section for the semi-inclusive DIS 



26 



bllowing 



27 



28|, 



da{ep-.^e'hX) ^ J d'k^d'p^^d'X^q{xs,k^:xsC)D{z,,p,^-X/ z,) 

xS{Xr, p)H{Q\xB, Zh] p)6^^\zhk± + pi^ + X± - p±) , (8) 

where q{xB, kj_), D{zh,pii_), S{\j_), and H are the transverse-momentum-dependent (TMD) 
quark distribution, fragmentation function, soft factor, and hard factor, respectively. We 
emphasize that the above factorization is valid in the leading power of p\/Q'^-, and all power 
corrections have been neglected.^ The energy dependent parameter C, C, and p have been 
introduced to regulate the light-cone divergences in the associated functions. In a special 
frame, we can simplify them as x^C^ = C'^ / zf^ = pQ^ [28]. The transverse momentum 
resummation can be performed by studying the evolution equation in terms of these vari- 
ables. We notice that, since the TMD quark distribution starts with nontrivial leading 
order expansion, the TMD fragmentation and soft factor in Eq. ([8]) are trivial at this order: 
D{zh,Pi±) = Dizh)5^'^\pi_L), S{\±) = S^'^\\±). However, at higher order, for example, the 
gluon radiation contribution to the SlDlS process in Fig. 1, we need to take into account 
nontrivial expansion of the fragmentation function and the soft factor up to order ^ 



The above factorization formalism was studied without considering the small-x resum- 



mation effects 



27 



281]. Here, we assume that the factorization argument can still hold 
when the hard momentum scale is much larger than the saturation scale and we 
can use the power counting method to study the leading contribution in this process. On 
the other hand, if Ql is the same order as (or even larger), the power counting used to 
argue the TMD factorization is no longer valid, and we will not have a TMD factorization. 
Similar studies for the heavy quark-antiquark production in pA {AA) collisions have also 



been discussed in [34]. 



^ The contribution from the TMD gluon distribution is power suppressed in term of p±/Q relative to that 
from the TMD quark distribution in Eq. ^ [3l and has been neglected in this equation. 
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FIG. 2: Transverse- momentum- dependent quark distribution calculated from small-x gluon splitting. 
The double line represents the gauge link contribution from the TMD quark distribution definition. 

As an important check, in the following we will compare the prediction from the TMD 
formula Eq. (8) to the dipole result Eq. (6) in the same kinematic region, ^ p'j_{Q'i)- To 
do that, we need to calculate the TMD quark distribution in nucleon/nucleus at small-x. 



This quark distribution is defined as 

72 



26| 



where P is the momentum for the hadron, x and k± are longitudinal momentum fraction of 
the hadron and transverse momentum carried by the quark. In the above equation^ C i s the 
gauge link introduced to guarantee the gauge invariance of the above definition 



26 



. At 



this particular order, the gluon splitting contribution to the TMD quark distribution can 
be calculated in the fc^-factorization approach at small x. We show the relevant Feynman 
diagrams in Fig. 2, where (b-d) diagrams come from the gauge link contributions. These 
diagrams have to be taken into account because the gauge field connecting to the hadron state 
(nucleon/nucleus) are dominated by the A'^ component in the fcrfactorization calculations. 
Their contributions are important to obtain a consistent and gauge-invariant result. The 
derivation is straightforward, and we have 

N f dr' f 

q{x, k^) = ^^j —j d%d'q^F{q^, x')A{q^, k^) , (10) 



where A has been defined in Eq. (7). This is the quark distribution calculated in the 
fct-factorization. In order to compare to the results we obtain above in the color-dipole 
formalism, we need to extrapolate in the leading logarithmic approximation at small-x, i.e., 
replacing the unintegrated gluon distribution F{qj_,x') by F{q±,x) in the above equation. 
Following this replacement, we will reproduce the differential cross section Eq. (6) calculated 
in the dipole framework at the leading order of p\/Q'^- Therefore, we have demonstrated 
that the small-x calculation of the differential cross section for the SIDIS process is consistent 
with the TMD factorization at this particular order. At even higher order, we will have to 
take into account the contributions from the fragmentation function and soft factor. At this 
order, they are trivial: D{zhiPii_) = D{zh)5^'^\pi±) and S{\j_) = 5^'^\\_l) where D{zh) is 
the integrated fragmentation function. We further argue that the TMD factorization will 
work at higher orders as well, because the power counting is valid when ^ p\{Ql) as we 
mentioned above. 

In the leading logarithmic approximation at small-x, we can further integrate out ^ in 
Eq. (10): 

fcx) = ^ / d%d-q^F{q^, X) {l - ^.^i^^'J^ly In J^^T^^ ' (^1) 



which is consistent with the result calculated before [23|]. A number of interesting features 



of this quark distribution have been discussed in the literature [23|, l2^. For example, in 
the small limit, the quark distribution saturates: xq{x^ki_)\k^^Q oc Nd'i'n'^; in the large 
k± limit, it has power-law behavior xg(x, oc Ql/k\. These two features will be 
manifest if we employ the GBW model for the unintegrated gluon distribution: -F(g_L,a;) = 
g-9i/Qa/7i-Qg, where is parameterized as = {x/ xoYGeV"^ with a;o = 3 ■ 10~^ and 
A = 0.28 (8|. Note that while the large behavior of the unintegrated gluon distribution 
F{q±, x) is incorrect in the GBW model (it falls exponentially instead of a power law), this 
bad feature does not translate to the TMD quark distribution: the convolution with the 
splitting kernel in Eq. (7) insures the proper leading- twist behavior. In Fig. 3 (left panel), 
we show the ratio of the TMD quark distribution xq{x, kx) relative to that at x = 10~^ as a 
function of k±_ for x = 10~^ and x = 10~^, respectively. From this figure, we can clearly see 
that the ratio remains unchanged when k± goes to 0, whereas the ratio is proportional to 
the ratio of at different x when k± is large. This clearly demonstrates that the transverse 
momentum dependence provides an important information on the saturation physics. We 
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FIG. 3: The transverse-momentum- dependent quark distributions at small-x: (left) at x = 10~^ 
and 10^'^ as ratios relative to that at x = 10~^ for fixed Q'^ = Qq = WGeV^ where the transverse 
momentum resummation effect is not important; (right) at different relative to that at = 
lOGeV^ for fixed x = 3 • 10"^ with Ql = IGeV'^ in the GBW model. 

have shown that these TMD quark distributions can be studied in semi-inclusive DIS process. 

Furthermore, the transverse momentum dependence is also sensitive to the QCD dynam- 
ics of the small-x evolution. In the above example, we took the simple parameterization 
from the GBW model This result shall be modified by the nonlinear evolution. For 
example, at large the unintegratedgluon distribution scales as {q^/Ql)'^" where Ac is 



the anomalous dimension 



18 



30 



35 



361]. In the DGLAP domain, we have Ac = 1 whereas 



in the BFKL domain it is Ac = 0.5. By solving the BK equation, it was found that Ac = 0.63 



for large rapidities F = In 1/x 30, 1371]. With this modification, the ratio of the TMD quark 
distribution at large will approach {QlY" instead of Ql- 

Another important QCD dynamics effect is the transverse momentum resummation 
which will affect the dependence of the spectrum. In the results we plotted in the left 
panel of Fig. 3, this effect was not considered, which correspond to the low results, for 
example, at = Qq = lOGeV^^. This effect can be studied by applying the CoUins-Soper- 
Sterman resummation method [s^]. There have been great applications of this method 



to various high-energy processes, in particular, in the semi- inclusive DIS at HERA 38 ] 
where important consequences have been observed. To demonstrate this effect in the TMD 
quark distribution at small-x we calculated above, we take the double leading logarithmic 
approximation (DLLA) to solve the evolution equation for the quark distributions. Under 
this approximation, we can write down quark distribution at higher in terms of that at 
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lower Ql [391, liOj^: 



(2^ 



0) ' 



(12) 



where the Sudakov form factor at the DLLA is defined by 



S{Q\Qlr) = \xi 



Ql 



An 



HQ'Qy) + Cor' 



(13) 



where we have also included a non-perturbative form factor contribution cor^lnQ^/Qo 41|. 
This resummation effect will shift the transverse momentum distribution to higher end when 

increases. As an example, in Fig. 3 (right panel) we show the typical changes for the quark 
distributions at Q"^ = 20, 50, lOOGeT^^ as compared to the lower Ql = lOGeV"^, with the 
following parameters: fixed coupling = 0.3, and Cq = 0.1 for the non-perturbative input 
for the form factor [ss!]. From this plot, we can see that indeed, the transverse momentum 
distribution becomes harder when is larger. We notice that neglecting the resummation 
effects in the left panel of Fig. 3 will introduce additional theoretical uncertainties in the 
predictions. However, we expect that much of the resummation effects will be cancelled out 
in the ratios of the quark distributions at the same Q^. 

In conclusion, we have studied the semi-inclusive DIS processes at small-x, and found 
that the quark distribution studied can be used as a probe of saturation physics. Especially, 
the ratio of the quark distributions is crucially depending on the saturation scale. We also 
studied the quark distribution at different values, and found that the resummation effects 
shift the distribution to larger k± with larger Q^. An ideal place to study this physics will 
be an electron-ion collider in the near future |42], where large nucleus target will provide 
an additional direction to study saturation. Meanwhile, we notice that the ratios plotted in 
Fig. 3 qualitatively agree with the experimental data from HERA 43|. Of course, in order 
to compare to these data, we have to take into account the fragmentation contributions 
to calculate the differential cross sections. We also notice that an extension to a study on 
the gluon transverse momentum distributions 4^ will have to consider both small-x and 
transverse momentum resummations. The result from this paper shall provide us confidence 
to carry out these important studies. 



^ Here, we approximate the energy dependent parameter ( in the TMD quark distribution by Q assuming 
they are of the same order 28 1. 



10 



F.Y. thanks Anna Stasto and George Sterman for interesting conversations a few years 
back about the topic studied in this paper. CM. would hke to thank the nuclear theory group 
at LBL for hospitality during the start of this work. We thank Al Mueller, Pavel Nadolsky, 
Jianwei Qiu, and Raju Venugopalan for stimulating discussions and comments. CM. is 
supported by the European Commission under the FP6 program, contract No. MOIF- 
CT-2006-039860. This work was supported in part by the U.S. Department of Energy 
under contract DE-AC02-05CH11231. We are grateful to RIKEN, Brookhaven National 
Laboratory and the U.S. Department of Energy (contract number DE-AC02-98CH10886) 
for providing the facilities essential for the completion of this work. 



[1] L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rept. 100, 1 (1983). 
[2] A. H. Mueller and J. w. Qiu, Nucl. Phys. B 268, 427 (1986). 

[3] L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994); Phys. Rev. D 49, 3352 

(1994); Phys. Rev. D 50, 2225 (1994). 
[4] see a review, E. lancu and R. Venugopalan, arXiv:hep-ph/0303204[ 



[5] I. Balitsky Nucl. Phys. B 463, 99 (1996). 

[6] Y. V. Kovchegov, Phys. Rev. D 60, 034008 (1999). 

[7] J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Nucl. Phys. B 504, 415 (1997); 
Phys. Rev. D 59, 014014 (1999); E. lancu, A. Leonidov and L. D. McLerran, Nucl. Phys. A 
692, 583 (2001). 

[8] K. J. Golec-Biernat and M. Wusthoff, Phys. Rev. D 59, 014017 (1999). 
[9] J. Bartels, K. J. Golec-Biernat and H. Kowalski, Phys. Rev. D 66, 014001 (2002). 
[10] E. lancu, K. Itakura and S. Munier, Phys. Lett. B 590, 199 (2004); G. Soyez, Phys. Lett. B 
655, 32 (2007). 

[11] K. J. Golec-Biernat and M. Wusthoff, Phys. Rev. D 60, 114023 (1999). 
[12] J. R. Forshaw, R. Sandapen and G. Shaw, Phys. Lett. B 594, 283 (2004). 
[13] K. J. Golec-Biernat and S. Sapeta, Phys. Rev. D 74, 054032 (2006). 
[14] C. Marquet, Phys. Rev. D 76, 094017 (2007). 

[15] H. Kowalski and D. Teaney, Phys. Rev. D 68, 114005 (2003); H. Kowalski, L. Motyka and 
G. Watt, Phys. Rev. D 74, 074016 (2006). 



11 



[16] J. R. Forshaw, R. Sandapen and G. Shaw, Phys. Rev. D 69, 094013 (2004). 

[17] C. Marquet, R. B. Peschanski and G. Soyez, Phys. Rev. D 76, 034011 (2007). 

[18] D. Kharzeev, Y. V. Kovchegov and K. Tuchin, Phys. Lett. B 599, 23 (2004). 

[19] A. Dumitru, A. Hayashigaki and J. Jahhan-Marian, Nucl. Phys. A 770, 57 (2006); A. Dumitru, 

A. Hayashigaki and J. JaliUan-Marian, Nucl. Phys. A 765, 464 (2006). 
[20] D. Boer, A. Utermann and E. Wessels, Phys. Rev. D 77, 054014 (2008). 
[21] J. JaUlian-Marian and Y. V. Kovchegov, Prog. Part. Nucl. Phys. 56, 104 (2006). 
[22] H. Weigert, Prog. Part. Nucl. Phys. 55, 461 (2005); D. N. Triantafyllopoulos, Acta Phys. 

Polon. B 36, 3593 (2005); G. Soyez, Acta Phys. Polon. B 37, 3477 (2006); S. Munier, Phys. 

Rept. 473, 1 (2009). 

[23] L. D. McLerran and R. Venugopalan, Phys. Rev. D 59, 094002 (1999); R. Venugopalan, Acta 

Phys. Polon. B 30, 3731 (1999). 
[24] A. H. Mueller, Nucl. Phys. B 558, 285 (1999). 

[25] The mean transverse momentum in this process was suggested to study the saturation physics, 
J. Qiu, Electron-Ion Collider Collaboration Meeting, Berkeley, December 11-13, 2008. 

[26] J. C. Collins and D. E. Soper, Nucl. Phys. B 193, 381 (1981) [Erratum-ibid. B 213, 545 
(1983)]; Nucl. Phys. B 194, 445 (1982). 

[27] J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B 250, 199 (1985). 

[28] X. Ji, J. Ma and F. Yuan, Phys. Rev. D 71, 034005 (2005). 

[29] A. M. Stasto, K. J. Golec-Biernat and J. Kwiecinski, Phys. Rev. Lett. 86, 596 (2001). 
[30] E. lancu, K. Itakura and L. McLerran, Nucl. Phys. A 708, 327 (2002). 
[31] C. Marquet and L. SchoefFel, Phys. Lett. B 639, 471 (2006). 

[32] F. Gehs, R. B. Peschanski, G. Soyez and L. Schoeffel, Phys. Lett. B 647, 376 (2007). 

[33] J. Kwiecinski, A. D. Martin and A. M. Stasto, Phys. Rev. D 56, 3991 (1997); K. Golec-Biernat 

and A. M. Stasto, larXiv:0905.1321l [hep-ph]. 
[34] F. Gehs and R. Venugopalan, Phys. Rev. D 69, 014019 (2004). 
[35] A. H. Mueller and D. N. Triantafyllopoulos, Nucl. Phys. B 640, 331 (2002). 
[36] J. L. Albacete, N. Armesto, A. Kovner, C. A. Salgado and U. A. Wiedemann, Phys. Rev. 

Lett. 92, 082001 (2004). 
[37] S. Munier and R. B. Peschanski, Phys. Rev. Lett. 91, 232001 (2003). 

[38] P. M. Nadolsky, D. R. Stump and C. P. Yuan, Phys. Rev. D 61, 014003 (2000) [Erratum-ibid. 

12 



D 64, 059903 (2001)]. 
[39] A. Idilbi, X. d. Ji, J. P. Ma and F. Yuan, Phys. Rev. D 70, 074021 (2004). 
[40] D. Boer, Nucl. Phys. B 806, 23 (2009). 

[41] G. P. Korchemsky and G. Sterman, Nucl. Phys. B 437, 415 (1995). 

[42] A. Deshpande, R. Milner, R. Venugopalan and W. Vogelsang, Ann. Rev. Nucl. Part. Sci. 55, 
165 (2005). 

[43] see, for example, C. Adloff et al. [HI Collaboration], Nucl. Phys. B 485, 3 (1997). 
[44] Y. V. Kovchegov and A. H. Mueller, Nucl. Phys. B 529, 451 (1998). 



13 



